The area generating function 
for simple- 2-column polyominoes 
with hexagonal cells 

Svjetlan Feretic * 
Faculty of Civil Engineering, University of Rijeka, 
Viktora Cara Emina 5, HR-51 000 Rijeka, Croatia 

Nenad Trinajstic ' 
The Rugjer Boskovic Institute, 
P. O. Box 180, HR-10 002 Zagreb, Croatia 

November 13, 2009 



Abstract 

Column-convex polygons were first counted by area several decades ago, and the 
result was found to be a simple, rational, generating function. In this chapter we 
generalize that result. Let a p-column polyomino be a polyomino whose columns can 
have 1, 2, . . . , p connected components. Then column-convex polygons are equivalent 
to 1-convex polyominoes. The area generating function of even the simplest gener- 
alization, namely to 2-column polyominoes, is unlikely to be solvable. We therefore 
define a class of polyominoes which interpolates between column-convex polygons and 
2-column polyominoes. We derive the area generating function of that class, using an 
extension of an existing algorithm. The growth constant of the new class is greater 
than the growth constant of column-convex polyominoes. A rather tight lower bound 
on the growth constant complements a compelling numerical analysis. 

1 Introduction 

The enumeration of polyominoes is a topic of great interest to chemists, physicists and 
combinatorialists alike 27 . In chemical terms, any polyomino (with hexagonal cells and 
with no internal holes) is a possible benzenoid hydrocarbon. In combinatorics, polyominoes 
are of interest in their own right because several polyomino models have mathematically 
appealing exact solutions. Furthermore, they are also relevant to various problems of tilings 
|25j . Polyominoes are extensively studied, in one form or another, in a wide variety of 
problems of interest to physicists. In particular, we note their investigation under the name 
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lattice animals, in the study of percolation [29, 33], in the graphical representation of the 
Ising model, and its extension to the Potts model, and in the study of the properties of 
branched polymers [2H1 1331 HE] • 

They are also a representative of a class of problems that appear to be unsolvable — 
notably the enumeration, by area or perimeter, of self-avoiding polygons, polyominoes and 
other classes of graphs that are relevant to the Ising model and to percolation. The principal 
line of attack on such problems is to simplify them until they are solvable, in the hope that 
the essential physics is not destroyed in the process. That is the approach taken in this 
chapter, where the model proposed, while still solvable, is closer to the ultimate problem of 
full polyomino enumeration than has previously been attained. Further, by use of Bousquet- 
Melou's [1] and Svrtan's [ST] upgraded version of the Temperley methodology [33], we give 
the solution of one problem previously out of reach due to its complexity. This development 
may spur further advances in this class of problem. 

One frequently cited polyomino model is that of column-convex polygons Q- There exist 
two main versions of column-convex polygons: the first composed of square cells and the 
second of hexagonal cells. Both versions have a rational area generating function. For the 
version with square cells, the area generating function was found independently by Polya [42) 
in 1938 or 1969 and by Temperley [44] in 1956. That was perhaps the earliest major result in 
polyomino enumeration. For the version with hexagonal cells, the area generating function 
was found by Klarner in 1967 [34] . The growth constant of square-celled column-convex 
polygons is /i = 3.205569 . . . , while the growth constant of hexagonal-celled column-convex 
polygons is /i = 3.863130 .... (By the growth constant we mean the limit linin^oo tfa^, 




where a n denotes the number of n-celled elements in a given set of polyominoes.) In both 
cases the area generating function is a simple pole, so that a n ~ const, x fi n . 

There exist some models which are supersets of column-convex polygons and are still 
solvable. Those models are called m-convex polygons [32], prudent polygons [23], cheesy poly- 
ominoes [IB], polyominoes with cheesy blocks [17], column- subconvex polyominoes [Tj5], and 
simple- 2- column polyominoes [18J . The former two models can be enumerated by perimeter 
and area, whereas the latter four models have been enumerated only by area. 

This chapter is a kind of companion to [18J. Namely, the enumeration of simple- 2-column 
polyominoes by area was first done in [18] , In the present chapter, we again enumerate 
simple-2-column polyominoes by area, but this time the cells are hexagons, whereas in [18] 
the cells were squares. As mentioned above, we make use of Bousquet-Melou's [3] and 
Svrtan's [21] upgraded version of Temperley's methodology [44]. The computations are 
rather long and intricate. Also, these computations are very similar to those which are 
done detailedly in [18 . Therefore in this chapter we only give an outline of the proof, 
though with enough detail that the method may be applied by others to new problems. 
Incidentally, in 18J, simple- 2-column polyominoes are called by their original name (coined 
by S. Feretic): simplex- duplex polyominoes. The name simple-2-column polyominoes was 
suggested subsequently (during the preparation of [5D]) by Tony Guttmann. 

In Section 2, we discuss the chemical relevance of polyominoes with hexagonal cells. In 
Section 3, we define the model. In Section 4 we give the formula for G(q,w), a generating 
function for simple-2-column polyominoes, in which the variable q is conjugate to the area 

1 We distinguish between polygons and polyominoes in that the former cannot have internal holes. As a 
consequence, the perimeter generating function for polygons has a non-zero radius of convergence, whereas 
for polyominoes the radius of convergence is zero. 
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and w is conjugate to the number of two-component columns of the polyomino. A truncated 
version of the proof is given in Section 5. In Section 6, we discuss the asymptotic behaviour 
of G(q, w), and give a tight lower bound on the growth constant, as well as a very accurate 
estimate. In Section 7 we conclude, outlining further work prompted by our results. 

Note that our solution essentially gives detailed information only about the area gen- 
erating function. The additional parameter w counts columns of a certain type. Thus we 
cannot give perimeter-area phase diagrams which are relevant to the description of vesicle 
collapse. Indeed, as is shown in [20], the perimeter generating function has zero radius of 
convergence (as is the usual case for polyominoes), which precludes such a phase diagram 
in its usual form. 

2 Chemical relevance of polyominoes 
with hexagonal cells 

Polyominoes with hexagonal cells can be used to model benzenoid hydrocarbons. Benzenoid 
hydrocarbons are a class of versatile conjugated organic molecules that have constantly been 
studied by a variety of researchers. This is so because benzenoid hydrocarbons found uses 
in experimental and theoretical research in chemistry 111! I26j , environmental chemistry 
[22l [40] , chemical technology [2j , cancer research [9[ [24j [31] , ecology [6] [37] , combinatorial 
chemistry [35] , graph theory (30] [45] , theory of aromaticity [39] , computational chemistry 
ESI SSI [47], polymer science [43], etc. 

Additionally, polyominoes and benzenoid structures under various disguises have been 
discussed regularly at the Dubrovnik MATH/CHEM/COMP Symposia. For example, this 
class of structures has been discussed there as fractal benzenoids [41] , as column-convex 
animals [TO] El [15] , in algebraic studies of benzenoids ,8 a , in leapfrog transformation [1] , 
maximum matchings and eigenvalues of benzenoid graphs 13J, in phototoxicity [12] . ring- 
currents in benzenoid hydrocarbons 38J, etc. 

Apparently research on various aspects of this class of structures continues to date and 
the present chapter is in this direction. However, before we proceed, we wish to point to 
perhaps the first scientific study on hexagonal structures. This study is due to Croatian 
polymath Rugjer Boskovic (1711-1787). He studied the structure of honeycomb. Hon- 
eycomb is made of hexagonal cells. Boskovic hypothesized that the hexagonal lattice is 
not chosen by chance. By using hexagonal lattice, bees economize on building material. 
Boskovic published this study, entitled De opium cellulis, in 1760. 

3 Definition of the model 

There are three regular tilings of the Euclidean plane, namely the triangular tiling, the 
square tiling, and the hexagonal tiling. We adopt the convention that every square or 
hexagonal tile has two horizontal edges. In a regular tiling, a tile is often referred to as a 
cell. A plane figure P is a polyomino if P is a union of finitely many cells and the interior 
of P is connected. Observe that, if a union of hexagonal cells is connected, then it possesses 
a connected interior as well, as a connected union of hexagonal tiles must be connected 
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Figure 1 : A column-convex polygon. 



through shared edges. Topologically, a connected union of square cells may be connected 
only at a shared vertex. Such unions are forbidden by the definition of polyominoes however. 

Let P and Q be two polyominoes. We consider P and Q to be equal if and only if there 
exists a translation / such that f(P) — Q- 

From now on, we concentrate on the hexagonal tiling. When we write "a polyomino" , 
we actually mean "a hexagonal-celled polyomino" . 

Given a polyomino P, it is useful to partition the cells of P according to their horizontal 
projection. Each block of that partition is a column of P. Note that a column of a polyomino 
is not necessarily a connected set. On the other hand, it may happen that every column of a 
polyomino P is a connected set. In this case, the polyomino P is a column- convex polygon. 
See Figure 1. 

By a 2- column polyomino, we mean a polyomino in which columns with three or more 
connected components are not allowed. Thus, each column of a 2-column polyomino has 
either one or two connected components. A simple- 2- column polyomino is such a 2-column 
polyomino in which consecutive two-component columns are not allowed. If c is a column 
of a simple-2-column polyomino, and c is a (left or right) neighbour of a two-component 
column, then c must be a one-component column. See Figure 2. 

4 The area generating function for 
simple-2-column polyominoes 
with hexagonal cells 

If a polyomino P is made up of n cells, we say that the area of P is n. Let 1Z denote the 
set of all simple-2-column polyominoes with hexagonal cells. 

In Theorem 2 below, we state a formula for the generating function 
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Figure 2: A simple- 2-column polyomino. 
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Theorem 1 The generating function G(q, w) is given by 



G(q,w) = 



NUM 
DEN' 



(1) 



where 



NUM = 



DEN 



(1 — q) A {a + 7 + 2arj — 2 7 e) + q 2 w(l — q) 2 {i + A — an — a/i 
+ (3l + (i\ - 7 k - 7y u + Si + S\ - 2e\ + 2ni + 2a(\ - 2ann 

- 2anu. + 2a9\ - 2/3e\ + 2[3ni + 2 7 e« + 2 7£y u - 2~/(l - 2j9l 

- 2Se\ + 25t]l) + Aq 2 w{l - q)(a\ - 7 t), 

(1 - qf{l - + S - e + n - a( + a9 + f3e - f3n + 7 C - 7<? 

- 5e + 8rj) - 2(1 - q) 3 (j + an - 7 e) 
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+ [ieu, - [3n^ - (39i + [39\ + 7O - 7 6>k - Sen + 5(l - 5(\ + Snn) 
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+ a(\ — otnn — anfi + a9\ — [3e\ + (inu + 7 eK + -feu, — j(i — j9l 

— SeX + 5ni) — 4q 2 w(ct\ — jt), 



{-iy^q 1 + 2i - 2 w % - 1 



£ 



(l- g )2i-2 . 


3) i - 1 q i2+3i w i 


4 2 
■(1-g*) 


(l-q) 2i - 

( 


"m=i(i-9 fc )" 

_g-ji-l^i 2 +2i-2 


4 q ' 
■{!-?) 

w % ~ 1 


(i-9) 2i - 2 -[ni=i(i-9 fe )" 


4 2 
■(1-qi) 



2i 



i—1 u 

q K 



2q* 



k=i 



Polyominoes with hexagonal cells 



7 



E 



(-3)" 



-l n i +3i w i 



^(i-5) 2i -[nt=i(i-9 fe ) 

q 3<f 



4 3 

■(!-«*) 



fe=i 



K = 



E 




^(i-?p- 2 -[nt=i(i-? fe )] -(i-^ 

1 + 4. ^ ^ 



2z 



+ 



fe=i 



-2* + l + 4-E7T 



„2fc 




■E 



(1 - 9 fc ) 2 (W) 2 _ 



(i-<z) 2j -[m=i(i 

2 *+ 4 -E 1 JL 1 3q 

k=i 

i-1 




q k )\ -(i -<?')' 

2 



^(l-5 fe ) 2 (W) 



In the above formulae, it will be noticed that (a) some of the numbers have an overline, 
and (b) no result is given for a, 5, e, 9, l, and fi. 

This is both to save space, and to highlight the close similarity between certain quantities. 
For all the quantities defined above, the overlines may be ignored. To obtain the formula 
for a from the formula for j3, replace 2 by 1. To obtain the formula for 6 from the formula 
for 7, replace the 3 by 4 and change (— 3) J_1 to (— 3) 1 . To obtain the formula for e from 
the formula for £, and also to obtain the formula for i from the formula for k, replace each 
of the 2's by 1. To obtain the formula for 9 from the formula for r], and also to obtain the 
formula for fi from the formula for A, replace each of the 3's by 4 and change (— 3) J_1 into 
(-3)*- 

5 Proof of Theorem 1 

Let P be a simple-2-column polyomino and let P have at least two one-component columns. 
Then we define the lower pivot cell of P to be the cell which is the lower right neighbour 
of the bottom cell of the second-last (i.e., second-rightmost) among the one-component 
columns of P. Wc also define the upper pivot cell of P to be the upper right neighbour of 
the top cell of the second-last among the one-component columns of P. 
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Let P be a simple-2-column polyomino and let P have at least one two-component 
column. Then we define the lower inner pivot cell of P to be the upper right neighbour of 
the top cell of the lower component of the last among the two-component columns of P. 
We also define the upper inner pivot cell of P to be the lower right neighbour of the bottom 
cell of the upper component of the last among the two-component columns of P. 

Observe that the lower pivot cell of a simple- 2-column polyomino P is not necessarily 
contained in P. The same holds for the upper pivot cell, the lower inner pivot cell and the 
upper inner pivot cell of P. 

Let S denote the set of those simple- 2-column polyominoes whose last (i.e., rightmost) 
column is a one-component column. It is convenient to first compute a generating function 
for the set S, and thence a generating function for all simple- 2-column polyominoes. So, let 



A(q,t,w) = J2 1^ ° f P ■ 1 C ~ 01 P ■ w 
Pes 



the height of the the number of 

component columns of P 



Next, we define three generating functions in two variables, q and w: Let A\ = A(q, 1, w), 
Bi = [Mteil and Ci = I . \ a2A ^y 

\ with t—1 L J with t—1 

The generating functions G(q, w) and A(q 1 i, w) are related by 

G(q,w)=A 1 + -^L_.C 1 . (2) 
Henceforth the notation A(q, t, w) will be abbreviated as A(t). 

In order to obtain a functional equation for the generating function A(t), we are going 
to suitably partition the set S. The blocks of the partition will be denoted S ai Sp, ... , S^, 
and the parts of A(t) coming from these blocks will be denoted A a (t), Ap(t), . . . , A^t), 
respectively. 

First, we define S a to be the set of those P <G S which have no other one-component 
column than the last column. We have A a (t) = j- 2 ^ + (j^rjjz^p ■ 

Let Sp be the set of those P e S \ S a which have the following two properties: the 
second-last column is a one-component column, and the last column contains the lower 
pivot cell of P. We have A (t) = 2 • A Y . 

Let S~{ be the set of those P e S\S a which have the following two properties: the 
second-last column is a one-component column, and the last column does not contain the 
lower pivot cell of P. We have A 7 (t) = ■ B x . 

Thus, SpUSy is the set of those P e S\S a whose second-last column is a one-component 
column. 

Let Ss be the set of those P e S \ S a which have the following three properties: the 
second-last column is a two-component column and the third-last column is (necessarily) 
a one-component column, the lower component of the second-last column and the third- 
last column have no edge in common, the lower pivot cell of P is contained in the upper 
component of the second-last column. We have Ag(t) = nr^q^TT^tp ' ^i- 

Let iS e be the set of P e S \ S a having the following three properties: the second-last 
column is a two-component column and the third-last column is a one-component column, 
the lower component of the second-last column and the third-last column have no edge in 
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common, and the lower pivot cell of P is contained in the hole of the second-last column. 
We have A e (t) = {1 - q % w . qt)i ■ A, (1 _ff^ gt)4 ■ A{qt). 

The definition of S^ is obtained from the definition of Ss by writing the word "upper" 
where the definition of Ss says "lower" , and by writing the word "lower" where the definition 
of Ss says "upper". We have A^(t) = Ag(t). 

The definition of S v is obtained when the changes just described are made to the defini- 
tion of S f (instead of to the definition of Ss). We have A v {t) — A f (t). 

Thus, Ss U S e U <J>£ U S v is the set of those P € S \ S a which, in addition to having a 
two-component second-last column, also have the property that P \ [the last column of P) 
is not a polyomino. 

Let Sg be the set of P e S \ S a which have the following three properties: the second- 
last column is a two-component column and the third-last column is a one-component 
column, each of the two components of the second-last column has at least one edge in 
common with the third-last column, and the last column also has at least one edge in 
common with each of the two components of the second- last column. We have Ag(t) = 

q 4 t 2 w p q 4 t 2 w a j q 4 t 2 w 4 ( ,\ 

(l-q)Hl-qt)' i ' 01 ~ (l-q)Hl-qt)± ' A ^ + (l-q)^(l-qtY ' > ' 

Let S L be the set of P e S \ S a which have the following four properties: 

• the second-last column is a two-component column and the third-last column is a 
one-component column, 

• each of the two components of the second-last column has at least one edge in common 
with the third-last column, 

• the last column has at least one edge in common with the lower component of the 
second-last column, but does not have any edges in common with the upper component 
of the second-last column, 

• the last column does not contain the lower inner pivot cell of P. 
WehaveA(t)= (1 _^_ gt) -C 1 . 

Let S K be the set of P e S \ S a having the following four properties: 

• the second-last column is a two-component column and the third-last column is a 
one-component column, 

• each of the two components of the second-last column has at least one edge in common 
with the third-last column, 

• the last column has at least one edge in common with the lower component of the 
second-last column, but does not have any edges in common with the upper component 
of the second-last column, 

• the last column contains the lower inner pivot cell of P. 
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A K (t) 



q tw 



q 4 t 2 w 



Ci 



(l-q) 2 (l-qtf (l- 9 )2(l-^)3 
q 4 t 2 w q 4 t 2 w 



' (1 - q) 2 (l qty ' Al (1 - q) 2 (l qt) 4 ' ^ ' 

The definition of 5a is obtained from the definition of S L by writing the word "upper" 
where the definition of S L says "lower" , and by writing the word "lower" where the definition 
of S L says "upper" . We have A\ (t) — A L (t) . 

The definition of 5 M is obtained when the changes just described are made to the defi- 
nition of S K (instead of to the definition of S L ). We have A^(t) = A K (t). 

Thus, Se U S L U S K U 5a U 5 m is the set of those P £ S\S a which, in addition to having a 
two-component second-last column, also have the property that P \ (the last column of P) 
is a polyomino. This means that S$ U S e U . . . U 6> M is the set of all P € S \ S a whose 
second-last column is a two-component column. 

The sets S a , Sp, . . . , 6> M form a partition of the set S. We have A(t) — A a (t) + Ap(t) + 
... + Afj,(t), and the expressions for A a (t), Ap(t),..., A^(t) are given above. Putting 
these things together, we get a functional equation for A(t). It is convenient to write that 
functional equation as 



A(t) = ^ 



qt 



1+Bt 



q 4 t 2 w 



2q 2 w 



•Ci 



qt 



(I -qt) Q 



{l-q) 2 {l-qtf 
3q 4 t 2 w 



■A x -B, 



A(qt). 



Ax 



2q 2 w 
3qH 2 w 



Ci 



(l-q) 2 (l-qty 



Ai 



(3) 



(l-q) 2 {l-qt) 4 

We solved equation ^ by iteration, as is usually done in the upgraded Temperley 
method. The iteration ended in 



(-3) i -y + 2i -H 2i - i w i - 1 



(_ 3 y-l q i>+2i-2 t 2 



+ E 



1 + Sl+ (T^F Cl 



i=l (i - g)ai-a . ^(1 - q H)\ ■ (1 - qH) 2 
(-3) i - 1 q e+3 H 2i w i 



E 



=i a-?) 21 - ni=i(i-^) -a-?**) 



( — 3) l q i +3i^2i,„i 



E- 



(4) 
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From equation 0$, we got a system of three linear equations in three unknowns: A±, B\ 
and C\. The first equation is just the case t = 1 of (0]). The second equation is obtained by 
differentiating ((4]) with respect to t and then setting t = 1. The third equation is obtained 
by differentiating (f4|) twice with respect to t, and then setting t = 1. 

Once the linear system is solved, relation ([2]) tells us how to obtain the sought-after 
generating function G. (To solve the linear system, we made use of the computer algebra 
package Maple.) 

5.1 A corollary to Theorem 1 

By setting w = 0, from Theorem 1 we obtain the well known result, discovered by Klarncr 

m- 

Corollary 1 The area generating function for column-convex polygons with hexagonal cells 
is given by 

q(l-q) 3 



G(q,0) 



1 - 6q + 10q 2 - 7q 3 



6 A bit of asymptotic analysis 

We write [q n ]f to denote the coefficient of q n in a power series / = f(q). 

Formula |T]) is very complicated, but still allows us to quickly compute Taylor polyno- 
mials of any reasonable degree. We actually computed the Taylor polynomial of degree 320 
for the function G(q, 1). It turned out that 



G(q,l) = g + 3<7 2 + ll( ? 3 + 44g 4 + 186g 5 + 806q 6 + 3518g 7 + 15349q 8 
+ 66797g 9 + 289960g 10 + 1256274q n + 5435860<? 12 + . . . . 

Then, for n — 2, 3,..., 320, we divided the coefficient [q n ]G(q, 1) by the coefficient 
[q n ~ 1 ]G(q, 1). This quotient gradually stabilizes, so that for n — 61, 62, . . . , 320 we have 



[g"]G(g,l) 
Vt-^Giq,!) 



= 4.322382971063 

rounded to 12 decimal places 



Simple concatenation arguments, first used by Klarner [34] . enable one to prove that 
H = lim I1 ^ 00 [<7 n ] 1 /™ = supjjg™] 1 /" . (The letter \x denotes the growth constant.) In this 
way, making use of [<7 320 ], we find that ijl > 4.294676, which is less than 1% below the best 
numerical estimate. Unfortunately, finding upper bounds is much more difficult. 

Next, for n = 1, 2, . . . , 320, we divided [q n ]G(q, 1) by 4.3223829710631654554™. Once 

again, the quotient stabilizes: for n = 54, 55, . . . , 320 we have 4 32238297i06 g 3i654554' 1 ~ 
0.127739087206. So, there is plenty of evidence that the coefficient [q n ]G(q, 1) (i.e., the num- 
ber of n-celled simple-2-column polyominoes) has the asymptotic behaviour [q n ]G(q, 1) ~ 
0.127739087206 • 4.322382971063™. The dominant singularity of G(q, 1) is a simple pole, 
located at 0.231353863527. 

For comparison, the area generating function for column-convex polyominoes is 
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G(q,0) = q + 3q 2 + llq 3 + 42g 4 + 162< 7 5 + 6 26 g 6 + 2419g 7 + 9346< 7 8 
+ 36106g 9 + 139483g 10 + 538841? 11 + 2081612g 12 + . . . . 

The number of n-celled column-convex polyominoes has the asymptotic behaviour 
[q n ]G(q,0) ~ 0.188419883819- 3.863130743243". The dominant singularity of G(q,0) is 
a simple pole, located at 0.258857405163 . 

As stated in [49], the area generating function for all polyominoes is 

q + iq 2 + llq 3 + Uq 4 + I86q 5 + 8Uq 6 + 3652q 7 + 16689g 8 
+ 77359g 9 + 362671g 10 + 1716033g n + 8182213g 12 + . . . 

and the growth constant of all polyominoes is 5.183147 .... Admittedly, the growth constant 
of all polyominoes is considerably greater than 4.322382 . . . , the growth constant of simplc- 
2-column polyominoes. 

7 Conclusion 

We have defined a class of polyominoes that interpolates between column-convex polyomi- 
noes and all polyominoes. The former have been solved, while the latter remain unsolved. 
The area generating function of column-convex polyominoes is just i q ]_ 7 i + i . On 

the contrary, the area generating function of our interpolating model (called simple-2-column 
polyominoes) is a cumbersome rational function of g-series. The area generating functions 
of column-convex polyominoes and of simple-2-column polyominoes both have a simple pole 
singularity, located at 0.258857. . . and 0.231353. . . , respectively. For all polyominoes, the 
corresponding singularity is at q = ^(polyomino) = 0.192932 . . . , and the singularity is of 
the form const, x |log(g c — q)\, rather than a simple pole [49] . We have also given the rigorous 
lower bound u > 4.294676, where [i is the growth constant of simple- 2-column polyominoes. 

In terms of possible extensions of this work, it is probably possible to compute the 
area generating function of simple- 2- columri2 polyominoes. Here, by a simple-2-column2 
polyomino we mean a 2-column polyomino in which runs of two consecutive two-component 
columns are allowed, but it is forbidden for three consecutive columns to each have two 
connected components. One reason for doing this is that the growth factor u is expected to 
be greater than that for simple-2-column polyominoes, and may set the benchmark in this 
regard. At present the situation is that for column-convex polygons the growth constant is 
u = 3.863130 . . . , for simple- 2-column polyominoes the growth constant is u = 4.322382 . . . , 
and for all polyominoes the growth constant is fi = 5.183147. . . . The polyomino model 
with a growth constant closest to the actual value for polyominoes is a directed model 
called multi- directed polyominoes [5] with a growth constant of fi sa 4.587894 . It would be 
interesting to compute the area of simple- 2-column 2 polyominoes to see if they had a growth 
constant closer still to that for polyominoes. 
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